Measurement of the form-factor ratios for D + —> K e + v e 
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Abstract 



We present a measurement of the form-factor ratios ry = V(0)/>li(0) and 
r 2 = A2(p)/Ai(0) for the decay D + — > K*°e + u e . The measurement is based 
on a signal of approximately 3000 D + — ► K*°e + v e , K *° — ► K~ir + decays re- 
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constructed in data from charm hadroproduction experiment E791 at Fermi- 
lab. The results are r v = 1.84 ± 0.11 ± 0.08 and r 2 = 0.71 ± 0.08 ± 0.09. 

PACS numbers: 13.20.Fc,14.40.Lb 
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Semileptonic charm decays are useful for probing the dynamics of hadronic currents be- 
cause Cabibbo-Kobayashi-Maskawa (CKM) weak mixing matrix elements for the charm sec- 
tor are well known from unitarity constraints. Form factors are Lorentz-invariant functions 
of q 2 , the square of the virtual W mass in the decay, that describe how strong interactions 
modify the underlying weak decay Form factor measurements in semileptonic decays 
test a variety of models and nonperturbative calculations. In addition, Heavy Quark Effec- 
tive Theory B relates form factors in semileptonic charm decays to those in bottom decays 
(at the same four- velocity transfer), which are needed to extract the weak mixing matrix 
elements \V u b\ and \V c b\ from semileptonic bottom decays. The vector form factor V(q 2 ) and 
the axial- vector form factors Ai(q 2 ) and A 2 (q 2 ) are relevant to the decay D + — > K*°e + i> e . 

Using data from charm hadropro duct ion experiment E791 at Fermilab, we reconstruct 
about 3000 D + — > K e + v e (and charge-conjugate) decays (three times the largest previous 
sample |J) and use the observed multidimensional distribution of kinematic variables to 
extract the form factor ratios r v = V(0)/Ai(0) and r 2 = A 2 (0)/A 1 (0). We assume the 
nearest-pole dominance model for the q 2 dependence of the form factors, F(q 2 ) = F(0) /(l — 
q 2 /m 2 ole ) where m po i e is the appropriate vector or axial-vector pole mass: my = 2.1 GeV/c 2 
or m A = 2.5 GeV/c 2 §. 

E791 is a hadroproduction experiment |J that generated charm using a 500 GeV/c 
tt~ beam incident on five thin targets (one platinum, four diamond) separated by gaps of 
about 13.6 mm. E791 ran with a loose transverse energy trigger and recorded 20 x 10 9 
interactions during the 1991-92 Fermilab fixed-target run. The important features of the 
E791 spectrometer for this analysis are the tracking system (23 planes of silicon microstrip 
detectors, 45 planes of drift and proportional wire chambers, and two large-aperture dipole 
magnets), two threshold Cerenkov counters that provide good K/tt separation over the 
momentum range 6-36 GeV/c, and a lead-liquid-scintillator electromagnetic calorimeter. 

In each event, we search for a candidate decay vertex (secondary vertex) with unit charge 
made from three charged tracks, separated from the reconstructed production vertex (pri- 
mary vertex) by at least 15a z where a z is the uncertainty in the longitudinal separation. 
The decay vertex is required to be at least one measurement error outside the nearest target 
material. One charged particle must be consistent with being an electron as determined 
by electromagnetic shower shape, the match between calorimeter energy and tracking mo- 
mentum, and the match between calorimeter and tracking position measurements. Electron 
identification efficiency is about 70%, while the probability for a pion to be misidentified as an 
electron is only 1 - 2%. One of the two remaining charged particles must have a Cerenkov 
kaon signature. Kaon identification efficiency is about 65% in the momentum range 6 - 
36 GeV/c and lower above this range. The probability for misidentifying a pion as a kaon 
is about 5% in the same momentum range and significantly lower above this range. Candi- 
dates consistent with being misidentified D + — > K~tt + tt + decays are removed. If electron 
and kaon candidates are oppositely charged, the decay is a candidate for D + — > K e + v e , 
K — ► K~n + . We call them "right-sign" decays. If electron and kaon candidates have 
the same charge, the decay is classified as "wrong-sign". The wrong-sign sample is used to 

model background in the right-sign sample. A clear excess of right-sign events compared to 

— *o 

wrong-sign events is seen in the Kn invariant mass distribution at the K mass. 

The final Kireve sample (see Fig. [I]) is optimized with a binary-decision-tree algorithm 
(CART H) that finds the set of "splits" in a set of single parameters or linear combinations 
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of parameters that best separates signal from background. We train CART using a sub- 
sample (~ 15%) of the right-sign candidates for signal and the wrong-sign candidates for 
background. CART selected a single cut involving a linear combination of four discrimina- 
tion variables: (a) separation significance of the candidate decay vertex from target material, 
(b) distance of closest approach of the candidate D momentum vector to the primary ver- 
tex allowing for the maximum kinematically-allowed miss distance due to the unobserved 
neutrino, (c) product over candidate D decay tracks of the distance of closest approach of 
the track to the secondary vertex divided by the distance of closest approach to the primary 
vertex, where each distance is measured in units of measurement errors, and (d) separation 
significance between the production and decay vertices. This selection criterion halved the 
number of wrong-sign events in the signal region, and kept 75% of right-sign events. 

Figure [I] shows mass distributions for the final right-sign (RS) and wrong-sign (WS) 
K7rev e candidates. The top left plot shows the distribution of M min , the minimum Kirev e 
mass kinematically allowed by the D direction as determined from the measured K, n, and e 
momenta and the positions of the primary and secondary vertices, with a 0.75 < M(Ktt) < 

1.05 GeV/c 2 cut for both the right-sign and wrong-sign events. The M m j„ distribution for 
true D + — ► K^Ti + e + v e events (with no detector smearing) would have a cusp at the D mass 
(1.869 GeV/c 2 ). The bottom left plot in Fig. |I] shows the M(Kn) distribution for events with 

1.6 < M min < 2.0 GeV/c 2 . Of these, there are 3595 right-sign and 602 wrong-sign events 
with 0.85 < M(K7r) < 0.94 GeV/c 2 (indicated by the vertical arrows in the figures), which 
are used to extract the form-factor ratios. The right-hand plot in Fig. [1] shows the difference 
between the right-sign and wrong-sign distributions. The net Kir signal is dominantly K*° 
as can be seen from the superposed fit of the M(Kir) spectrum to a pure Breit-Wigner shape 
with the mass and width fixed to the known values of the K*° resonance. There is an excess 
of events over that expected for a Breit-Wigner distribution in the region below the M(Ktt) 
range used in the analysis. The assumption that wrong-sign events accurately model the 
size and shape of the right-sign background is addressed in the discussion of systematic 
uncertainties below. 

The kinematic variables used in extracting the form factor ratios are the square of the 
invariant mass of the virtual W (q 2 ) and three angles. The polar angle 8 e , measured in 
the virtual W (or ev) rest frame, is the angle between the charged lepton and the direction 
opposite the K*°. The polar angle 9y, measured in the K*° rest frame, is the angle between 
the K and the direction opposite the virtual W . The azimuthal angle x is the angle between 
the momentum projections of the electron and K in the plane perpendicular to the K*° 
direction in the D rest frame. To calculate these variables, the neutrino momentum is 
estimated up to a two-fold ambiguity from the D flight direction as determined by the 
measured positions of the D production and decay points, and the measured momenta of 
the charged decay products. Monte-Carlo simulation shows less kinematic variable smearing 
for the solution which results in the lower D momentum, so it is used. 

We extract the form factors using an unbinned maximum-likelihood method that uses a 
Monte-Carlo simulation in the evaluation of the likelihood function The production 

physics and detector response are simulated for an event set that is passed through the same 
analysis chain as the data. The simulated events are generated with specified form-factor 
ratios (0.82 for r 2 , 2.00 for ry) 0. The likelihood of the data sample is calculated, for 
any given set of theoretical parameters, by computing the density of Monte-Carlo events 
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in a specified volume around each data point, where the simulated events are distributed 
according to the theoretical parameters under consideration. To avoid generating separate 
Monte-Carlo samples for every set of theoretical parameters considered in the fit, a single 
Monte-Carlo sample is reweighted so that the weighted events give the correct density about 
each data point. As long as the Monte Carlo accurately simulates both the charm produc- 
tion process and the detector response, acceptance and smearing effects are automatically 
incorporated into the fit. The wrong sign candidates are used to incorporate backgrounds 
into the fit with a similar technique. We developed a second method to extract form factors 
that keeps both neutrino momentum solutions and we use it as a check. To account for the 
wrong solution, we use Monte Carlo simulation to determine a feedthrough matrix that gives 
the probability that an event appearing in one region of the space of measured kinematic 
variables would appear in another region when the other solution for neutrino momentum is 
used. This matrix and the observed distribution of data events (both solutions) determine 
the fraction of data events that correspond to the correct neutrino-momentum solution in 
each region of kinematic-variable space. Each fraction is then used in a binned maximum 
likelihood fit, with background modeled as in the first method. 

The first fitting technique gives r v = 1.84 ± 0.11 ± 0.08 and r 2 = 0.71 ± 0.08 ± 0.09 with 
a correlation coefficient of —0.13. The first error is statistical, and the second systematic 
(discussed below). The above results have been corrected for biases of less than 10% due to 
the technique, determined from fitting Monte-Carlo samples with known form-factor ratios. 
The second fitting technique gives ry = 1.78 ± 0.12 and r 2 = 0.68 ± 0.07, where the errors 
are statistical only. The Monte Carlo indicates no correction for bias is required. The results 
from the second fitting technique are consistent with the first. 

Possible sources of systematic uncertainties were considered and the most important are 
summarized in Table |I[ To estimate the effects of possible inaccuracies in the Monte-Carlo 
simulation of the detector response, 15 different sets of selection criteria were generated 
using different training samples for CART. The spread in the resulting form-factor ratios 
gives our estimate of this systematic uncertainty [HJ. We estimate the uncertainty due 
to our modeling of the background by varying our assumptions about the amount and 
the distribution of the background in the four-dimensional kinematic variable space. We 
have determined that D* + — > D°7i + ,D° — > K~e + v e , and D + — » K* Q n°e + v e modes do not 
contribute significantly to the background in the signal region. Other sources of systematic 
uncertainty are the limited size of the Monte-Carlo sample and uncertainties in charged- 
particle identification efficiency. The contributions from each source are added in quadrature 
giving the total error of 0.08 and 0.09 for ry and r 2 , respectively. 

Figure ^| shows the projections of the kinematic variables cos 9y, cos9 e , and x f° r data 
and for Monte-Carlo events that have been weighted according to the best fit values for the 
form-factor ratios. To reveal the correlations between the observed kinematic variables, we 
show plots for each variable for two ranges of a second variable. The confidence level for the 
consistency of the Monte Carlo and data projections is shown on each plot. 

Throughout, we have assumed the nearest pole dominance q 2 dependence F(q 2 ) = 
F(0)/ (l—q 2 /m 2 A v ) which can be approximated by a form linear in q 2 , F(q 2 ) = F(0)(l+p 2 q 2 ). 
This is a valid approximation in the accessible q 2 range (0 - 0.947 GeV 2 /c 4 ). We per- 
form a three-parameter fit for the slope p 2 A and the form-factor ratios ry and r 2 (fixing 
Pa/ Pv = m y/ m i)- The result for r v is 1.88 ± 0.11, which is close to the two-parameter 
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fit. The results for r 2 and p\ are strongly anticorrelated, so the statistical uncertainty on 
these two parameters is large. The result for r 2 is 0.98ig;}5, which is about two standard 
deviations higher than the result of the two-parameter fit. The result for p\ is — 0.06ig;og 
GeV _2 c 4 , which is about two standard deviations lower than the theoretical expectation 
p\ = l/ m \ = 0.16 GeV~ 2 c 4 . 

Table |I| compares the form-factor ratios measured by this experiment and previ- 
ously published results. All measurements are in accord. Using the world averages for 
B(D + — > K e + v e ) and the D + lifetime ||, we extract the values of the form factors A%, A%, 
and V at q 2 = and at q 2 = q 2 max . We account for both the finite width of the K 
and the correlation between the measured form-factor ratios |[12|| . Table [TTI| compares these 
results with predictions from Lattice Gauge calculations Ji3Hl6| and a quark model calcu- 
lation jnj based on Heavy Quark Effective Theory. The former are in fair agreement with 
the experimental results, while the latter is not, in particular for A 2 (q^ nax ). 

In summary, we have used a sample of 3000 signal events to extract the form-factor ratios 
in the decay D + -> K*°e+u e : r v = 1.84 ± 0.11 ± 0.08 and r 2 = 0.71 ± 0.08 ± 0.09. The 
combined statistical and systematic uncertainties are a bit less than half the best previous 
measurement. The form- factor ratios are important for improving our understanding of 
the dynamics of hadronic currents and might improve our ability to extract CKM matrix 
elements from B semileptonic decays. 
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TABLES 



TABLE I. Contributions to the systematic uncertainty. 


Source of Uncertainty 


u rv 




Simulation of detector effects 


0.03 


0.03 


Monte Carlo volume size 


0.04 


0.04 


Background volume size 


0.04 


0.05 


Amount of background 


0.02 


0.05 


Particle identification efficiency 


0.05 


0.01 


Fitting technique 


0.01 


0.01 



Total Estimate 0.08 0.09 



TABLE II. Results of this analysis (E791) and comparison with other experiments. The ap- 
proximate number of signal events and the lepton type used in each analysis is also listed. 



Experiment 



r v = V(0)/A 1 (0) 



r 2 =A 2 (0)/A 1 (0) 



Events 



E791 
E687 @ 



E653 [DJ 
E691 § 



1.84 ±0.11 ±0.08 
1.74 ±0.27 ±0.25 
2.00±g;|f±0.16 
2.0 ±0.6 ±0.3 



0.71 ± 0.08 ± 0.09 
0.78 ±0.18 ±0.10 
0.82±g;|§ ±0.11 
0.0 ±0.5 ±0.2 



3000 (e) 
900 O) 
300 fa) 
200 (e) 



TABLE III. Form factors extracted from the measured ratios ry and r 2 , and several theoretical 
predictions. 



Group 



Ai(0) 



A 2 (0) 



V(0) 



E791 
APE 
Wuppertal |L| 
UKQCD 



ELC g6| 



0.58 ±0.03 
0.67 ±0.11 



0.61 



+0.11 
-0.09 



o 7n+ - 07 

u - /u -0.10 

0.64 ±0.16 



0.41 ±0.06 
0.49 ± 0.34 



0.83 



+0.23 
-0.22 



0.66l°;{° 
0.41 ±0.28 



1.06 ±0.09 
1.08 ±0.22 

1 Q4+0.31 

0.86 ±0.24 



^1 (.Qmax) 



E791 

ISGW2 p| 



0.68 ±0.04 
0.70 



0.48 ± 0.08 
0.94 



1.35 ±0.12 
1.52 



S 



FIGURES 




FIG. 1. Final sample of D + — > K e + u e , K* — > K tt + candidates. The distributions are 
described in the text. 



(a) (b) 




COS^v COST? e X 

FIG. 2. Data distributions (black dots with errors) overlaid with Monte Carlo (histogram) for 
(a) cos^y for q 2 /q 2 nax ^ 0-5 (top) and q 2 /q 2 nax > 0-5 (bottom), (b) cos# e for q 2 /q 2 nax < 0.5 (top) 
and q 2 /q^nax > 0.5 (bottom), (c) % for cos^y < (top) and cos6V > (bottom). The Monte-Carlo 
events are weighted according to the best-fit values of the form-factor ratios. 
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